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ABSTRACT 

We derive expressions for the total and linearly-polarized synchrotron emissivity of an 
element of plasma containing relativistic particles and disordered magnetic field which 
has been sheared or compressed along three inde pendent directions. Our treatment 
follows that given by Matthews & Schcuer (1990) in the special case of a power-law 
electron energy spectrum. We show that the emissivity integrals depend on a single 
parameter, making it straightforward to generate one-dimensional look-up tables. We 
also demonstrate that our formulae give identical results to those in the literature in 
special cases. 
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1 INTRODUCTION 

The inference of the three-dimensional structure of the mag- 
netic field in a source of synchrotron radiation from the 
observed emission at a single frequency is not unique, but 
considerable progress can be made if well-resolved images 
of total intensity and linear polarization are available. We 
have developed models of jets in low-luminosity radio galax- 
ies which pre dict the observed e mission for comparison with 
observations ( Laing et al. 199S| ; Laing & Bridle, in prepara- 
tion) . Our fundamental assumption about the field structure 
is that it is random on small scales, but made anisotropic by 
shear and compression. Given the large number of observed 
data points, the need to integrate along the line of sight 
and the iteration required to optimize the fit between model 
and data, we found that existing algorithms for the calcu- 
lation of synchrotron emissivity from disordered fields were 
too cumbersome. In order to minimize the resulting compu- 
tation, we have developed an efficient method of evaluating 
the Stokes parameters I, Q and U for synchrotron radiation 
from electrons radiating in a three-dimensional, anisotropic 
magnetic field. This technique should be applicable to a va- 
riety of astrophysical synchrotron sources and is the subject 
of the present paper. 

The basic theory of synchrotron emission from an 
optically-thin power-law energy distribution of relativistic 
electrons n(E)dE oc E~^^°''^^^ dE is well-docmn ented (e.g. 
Rybicki fc Lightman 1979[ [Hughes fc Miller 199l| ). The spec- 
trum of the emitted radiation also has a power-law form 
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with spectral index a {S oc p""). For a uniform field, 
B, at an angle 9 to the line of sight, the emissivity is 
(X [B sinOY^" and the degree of polarization has its maxi- 
mum value po = (3Qf-l-3)/(3a-l-5). The observed E- vector is 
perpendicular to the projected field direction on t he sky in 
the absence of propagation effects. Scheuer (1967) pointed 
out that the same results hold for a one-dimensional disor- 
dered field with many reversals: /, Q and U depend only on 
the orientation of the field, not on its vector direction. 

Expressions for synchrotron emission from a two- 
dimensional disordered field with equal component variances 
along two orthogona l directions and no component in the 
third were d erived by Laing (1980), and discussed further by 
Laing (1981). The more general case of an initially isotropic 
field compressed by an arbitrary fact or along one direction 
was considered by [Hughes et al. (19*85 ) . These authors wrote 
down expressions for elements of field in terms of spherical 
polar angles and averaged over the angles to derive the emit- 
ted Stokes parameters. For a field structure described by one 
compression factor, this approach is well suited to numerical 
calculation, but it be comes extremely cumbersom e if veloc- 
ity shear is included. Matthews & Scheuer (1990, hereafter 



MS) showed that a much simpler treatment is possible if it 
is assumed that the initial field is isotropic with a Gaussian 
distribution of field strengths. Their formulae also include 
the effects of inverse Compton and synchrotron losses, which 
cause a deviation of the emitted spectrum from a power-law 
form. In this paper, we describe a simplified version of their 
method, which applies to the three-dimensional case, but is 
restricted to power-law spectra. The advantage is that com- 
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putation of the emissivity is reduced to a one-dimensional 
integration or table look-up. 

The new formulae are derived in Section ^ In Section ^ 
we illustrate the use of the method in practice by rederiving 
expressions in the literature for a number of special field 
configurations. Section briefly summarizes our results. 



COMPUTATION OF THE STOKES 
PARAMETERS 



+ 2{ax ay + bxby + cxCy) sin xo cos xo] (3) 

= Bo[{ax + b% + c%) sin^ Xo 

+ {ay + by + cy) cos^ Xo 

— 2{axaY + bxby + cxcy) sin xo COS xo] (4) 



Here, Bo is the rms field along any direction for the initial 
isotropic field. Note that MS, in their equation (3.8), use the 
mean square value of the total field {3Bq in our notation). 



2.1 Assumptions 

The jet emission is taken to be optically-thin synchrotron 
radiation with a power law frequency spectrum of index 
a as defined earlier. We make the conventional assump- 
tions that the pitch-angle distribution of radiating electrons 
is isotropic. The observing frequency is taken to be high 
enough that Faraday rotation may be neglected, or corrected 
accurately. 

The magnetic field is assumed to be disordered on small 
scales with zero mean. We and others have argued that this 
is likely to be a good approximation for extragalactic radio 
sources, except perhaps on the smallest lin ear scales (Laing 
1981; Begelman, Blandford & Rees 1984; Laing & Bridle, 



in preparation). We note, however, that if one of the three 
field components is vector-ordered, then the results of our 
emissivity calculations are unaltered. 



2.2 Geometry 

Following MS, the field structure is assumed to be generated 
from an element of plasma containing radiating particles and 
an initially isotropic field distribution by shear and compres- 
sion. We start with a small cube of plasma whose sides are 
defined by three unit vectors in a coordinate system x, y, z 
moving with the emitting plasma (the distinction between 
emitted and observed frames is important for relativistic 
flow). The coordinate axes are usually chosen along obvi- 
ous natural directions. This element of particles and field 
is then compressed and sheared, forming a parallelepiped 
whose sides are defined by the vectors a, b and c. We then 
calculate the components of these vectors in a coordinate 
system X, Y, Z fixed in space, with Z along the line of sight. 

MS show that the position angle xo of the appar- 
ent magnetic field (perpendicular to the E-vector for syn- 
chrotron radiation) is given by: 



tan 2x0 



2{axay + bxby + cxCy) 



+ b\ 



for any spectral index, ax is defined by: 
ax ~ a.X/a.b x c 



(1) 



(2) 



where X is a unit vector along the X-axis, and so on. We 
will refer to ax ,aY ,bx ,by , cx ,cy collectively as the edge 
vector components. We rotate about the line of sight by an 
angle xo into a coordinate system X, Y , Z {Z — Z) in which 
Stokes U — 0. The rms field components along the X and 
Y axes are Bj^ and By, where: 

B% = Bo[{ax + bx + Cx) cos'^ Xo 
+ {ay + by + Cy) sin^ Xo 



2.3 Gaussian field distributions 

For a power-law electron energy distribution, the observed 
emissivity from a element of field of strength B is always of 
the form B^'^'' x a geometrical factor which is independent 
of B. The simplest initial field distribution is one where the 
field has constant intensity and random direction. Whilst 
this gives simple expressions for the emissivity if a = 1, the 
more general case is extremely messy. MS showed that a 
much simpler treatment is possible if the initial distribution 
of field strengths is isotropic and Gaussian, with probability 
distribution function (PDF): 

F{B^,By,B,)dB^dBydB, 



= (27rBo)~^''^exp 



B^ 
2Bi 



dB^dBydB^ 



(5) 



where Bo is again the rms field in one coordinate, Bx, By and 
Bz are the fields along x, y and z and B^ — B^ + By + B^. 

In this case, the PDF of the field after shear or com- 
pression can still be described by a multivariate Gaussian 
function. The PDF of the components of magnetic field per- 
pendicular to the line of sight (responsible for synchrotron 
emission) then has its major and minor axes along X and Y, 
with the rms fields Bj^ and By given earlier. Unless By ~ 0, 
the joint PDF for B± and the position angle x is as given 
by MS: 



Fx±{Bi 



,x)dB±dx 
B_ 



2TvB^By 



exp(— ui?x — vB^ cos2x)dBj^dx 



where 



(6) 

(7) 
(8) 



The PDF for B± is given by integration of equation (^) over 
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f±(Bx) 



Bx exp(— 
2nB^By 

B± exp{-uBi)Io{-vBl 
BxBy 



exp(— uBx cos 2x)dx 



(9) 



In is the modified Bessel function of or der n and we have 
used the in tegral representation given by Abramovitz fc Ste- 
gun (197c| ), equation (9.6.16). This is identical to equation 
(3.12) of MS, since Jo is an even function. 
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2.4 Synchrotron emissivities for a Gaussian field 
distribution 

We work with total and polarized intensity functions ji and 
jp defined such that the rest-frame emissivities are given 
by ej = k{n,a)ji and ep = k{n,a)jp where k{n,a) is pro- 
portional to the particle density n and the functions are 
normalized to have unit value for a uniform field of unit 
strength perpendicular to the line of sight, ji and jp are 
then functions of the field strength, geometry and spectral 
index: 



Ji = 



Jp 



Bx^Y Jo 



B^+°'exp{-uBl)Ioi-vBl)dB^ (10) 



'0 

Po 



po / dSx / dxiJx "cos2xi^xx(Sx,x) 



B^^" exp{-uB'i)Ii{-vBl)dB^ (11) 



W e have again used the integra l representation for Ii given 
by Abramovitz fc Stegun (197C| ), equation (9.6.19). It is con- 
venient to write these formulae in terms of the dimensionless 
variables q = B± /Bj^ and rj = B^ /Bj^ (the axial ratio of the 
field PDF) and to define the functions Go{q) = exp{—q)Io{q) 
and Gi(g) = exp(— g)/i(g). For q :> n, 



In{q) 



: exp(g) 



so G„{q) oc q The total and polarized intensity func- 
tions then become: 



JI = 



JP 



■exp(-gV2)Go 



^ ^-1 



exp( 



-gV2)Gi 



^2 



dq (12) 



dq (13) 



In the special case By =0 (?7 = 0), the PDF for B± is a 
one-dimensional Gaussian function: 



F^{Bi_)dBi_ = 



.Bl 



1/2 



exp 



Bi_ 
'2BI 



dB, 



and the integral for ji can be solved explicitly: 
JI = li-^ -n -^2^^+"^'^r(l + a/2) 



JP 



^l + a^-l/2r,(l + a)/2-f 

Poji 



(14) 



(15) 
(16) 



using the relations Iq{w) — Jo(iw) and Ii{w) 
and the substitution w = q^ . The results are: 

ji 



jp = poB^{l~Tl 



-iJi (iw) 

(17) 
(18) 



and remain valid if = 0. The degree of polarization is 
simply: 



P = jp/ji = Po 



1-rf 
1 + 



(19) 



Finally, the Stokes parameters in the fixed X, Y coor- 
dinate system are given by: 



JQ = jp cos 2x0 
ju = jp sin 2x0 



(20) 
(21) 



Note that these definitions of Q and U refer to the direction 
of the apparent magnetic field, and that the position angle 
of the synchrotron E- vector is xo + n/2. 



2.5 Emissivities for an initial field of constant 
magnitude 

An initially isotropic Gaussian field with component rms 
-Bo can be regarded as the superposition of a set of isotropic 
fields of constant magnitude B with a distribution: 



Fw{B)dB 



(27r)i/2B3 



B^ exp( 



-B^/2Bl)dB 



(22) 



Consider two cubes, both containing initially isotropic fields, 
but with different distributions. The first has a Gaussian 
distribution, with component rms Bq, as before. The second 
has a field of constant strength Bo. They are stretched and 
sheared in identical fashions. Since the geometrical factors in 
the emissivity calculation are independent of field strength, 
the ratio of emissivities (Gaussian/constant) is: 



B'+°'Fw{B)dB 



-°o 



-,(5+c.)/2 -1/2 



exp( 



^t^^r'+'^dt 



2(3+a)/2^-l/2p(.2 + a/2) 



(23) 



again using GR 8.310. 

The total and polarized emissivities for the case of a 
constant initial field are ji /r and jp /r, respectively, r = 3 if 
Q = 1, refiecting the facts that the emissivity is proportional 
to the square of the field and that the mean square fields 
along a given direction are Bq and -Bo/3 for the Gaussian 
and constant-field cases, respectively. 



using the integral representation of the Gam ma func- 
tion, r(w), given in Gradshteyn & Ryzhik (1994, hereafter 



GR) 8.310. In this situation, there is no field component 
along the Y axis and the degree of polarization takes its 
maximum value. 

If a = 1, there are simple analytical solutions. The inte- 
grals for 77 7^ can be written in the form given by GR 6.623 



2.6 Implementation 

The integrals can be evaluated usi ng standard rout ines, for 



example QROMO and midexp from Press et al. (1992 ). These 
implement Romberg integration for an integrand which de- 
creases exponentially and an infinite upper integration limit. 
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The functions Go and Gi can be computed by slightly modi- 
fied versions of th e subroutines bessiO and BESSll, also from 
Press et al. (1992). The emissivity functions depend only on 
the ratio of the field rms's, rj, for a given spectral index. 
Look-up tables can be generated as functions of rj at the 
beginning of a computation and thereafter polynomial in- 
terpolation can be used to derive values as required. 

The emissivity functions ji and jp are plotted against 
axial ratio rj in Fig. Q for a range of spectral indices typical 
of optically-thin synchrotron radiation. The corresponding 
degrees of polarization, p — jp/ji, are shown in Fig. ^ Al- 
though the variation of total intensity with spectral index is 
substantial, that of the degree of polarization is not, and the 
analytical formula for a — 1 (equation |l9| ) provides a good 
approximation f or most sour ces, a s previously noted in less 
general cases by Laing (198C) and Hughes et al. (1985). 



3 EXAMPLE CALCULATIONS 

In this section, we outline the use of the method in practice 
and demonstrate that the formulae given in the previous 
section are equivalent to expressions in the literature for a 
number of simple field distributions. These distributions are 
non-Gaussian, and their emissivity functions are denoted by 
// and fp, to distinguish them from the Gaussian equiva- 
lents. 



3.1 General approach 

In one class of problem, we are given the rms field compo- 
nents along three orthogonal directions. Such a configura- 
tion can always be generated from a unit cube containing 
an isotropic field by stretch or compression along the three 
axes. We can choose orthogonal vectors a, b and c along the 
axes so that the volume V of the resulting cuboid has any 
specified value. The initial field Bo and the magnitudes of 
the vectors are determined by: 

{Bl}^^^ = aBo/V (24) 
bBo/V (25) 
cBo/V (26) 
abc (27) 



V 



so we can solve for Bo, a, h and c. 

We can extend this approach to follow the emission from 
an element of fluid in a prescribed flow pattern (including 
the effects of shear) by tracking the vectors a, b and c as in 
MS, using the values given above as initial conditions. 



3.2 One-dimensional fields 

To get a one- dimensional field, we start with a cube of ini- 
tially isotropic, Gaussian field of rms Bo in each coordinate 
and expand it by a factor c — > oo in one direction (z), leav- 
ing the other two dimensions constant. The rms field in the 
z direction remains Bo, whilst those in the other two co- 
ordinates — > 0. Suppose that the 2-axis makes an angle 6 
with the line of sight. In the formulation of Section]^, 77 = 
and the only non-zero edge vector component is cx = sin 6. 
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Figure 1. Plots of the emissivity functions for total intensity 
(lower panel) and linear polarization (upper panel) against the 
axial ratio of the magnetic field PDF, r\ for different spectral 
indices a. The values of ct are: 0.0 (dotted), 0.5 (short dashes), 
1.0 (full line), 1.5 (long dashes) and 2.0 (dash-dot). 



The emissivities are given by equations (|l^) and ( [l^ ) with 
Bj^ = Bo sin 0. 

Alternatively, we can use the standard emissivity func- 
tions for a one- dimensional field of strength B: 



fi 
fp 



{Bsme)^+" 
Pofi 



(28) 
(29) 



Synchrotron emission 5 



0.8 




0.2 0.4 0.6 0.8 
Axial ratio 77 



Figure 2. Plot of the degree of polarization against the axial 
ratio of the magnetic field PDF, rj for different spectral indices 
a. The values of a are: 0.0 (dotted), 0.5 (short dashes), 1.0 (full 
line), 1.5 (long dashes) and 2.0 (dash-dot). 



together with the PDF for the total field, which is a one- 
dimensional Gaussian function: 

1/2 



Fm{B)dB = 



exp 



'2Bl 



dB 



(30) 



Integration of B^^" over the Gaussian PDF gives equa- 
tions (|l|) and (^ with B^ = Bo sin 

3.3 Two-dimensional fields 

Analytical expressions for the total and polarized emissiv- 
ities for a sheet of disordered field with no component in 
one direction but equ al component s in the other two coor- 
dinates were given by Laing (198C), who considered a field 
of constant strength B and random direction, confined to a 
plane (note that this distribution is not produced by simply 
squashing a three-dimensional field of constant magnitude 
and random direction). 

In order to compare these formulae with their equiva- 
lents from Section ^, we consider a cube of initially isotropic, 
Gaussian field of rms Bo in each coordinate which is ex- 
panded by factors a, a and 1/a along the x, y and z axes in 
the Umit a ^ oo, so that there is no field in the z direction, 
but the rms values in the x and y directions remain equal to 
Bq. Since the z-axis makes an angle 6 with the line of sight, 
the edge vector components are: 



ax = 1, fox = 0, cx = 
ay = 0, foy = cos 6, cy =Q 



(31) 



given by inserting these values into equations ( p^ ) , ( p^ ) , ( p^ 
and (nq). For 6 / 7r/2, we make the substitutions: 



1 + cos^ e 



2 cose 



(32) 
(33) 



into the expressions for ji and jp , so 

Dl+a„2+<:«C0S^^"6' 



sm 



oo 



l + o 



X / t~ exp[-t™(?j;^ - l)~''^]h{t)dt (34) 



2 + a a 



Dl + «r,2-|-Q, COS 

JP = poBo 2 



sm 



e^p[-tw{w'^ ~ l)'^'^]h{t)dt (35) 



The integrals can be written in terms of Associated Legendre 
and Gamma functions using GR 6.624 and 8.731: 

3. = Bo^-2H-r(^)(cose)^ 



X Pc + l 



1 + cos^ e 

2 cos 6* 



jp = poBo^-^"2 — r 



X Pl + 1 



1 + COS'' I 

2 cos 6* 



(36) 



(37) 



unless 6 — 7v/2, in which case equations (|l^) and ( p^ apply, 
with Bj^ = Bo. 

For a spectral index a — 1, equations (^^ and ( p^ 



ji = Bo(l + cos^( 
jp = poBq sin^ e 



(38) 
(39) 



Alternativ ely, we can start from the expressions given 
by Laing (1980) and average over the PDF for the total field. 
We make three modifications to the notation: 

(i) We normalize over the PDF for azimuthal angle, in- 
troducing an extra factor of 27r. 

(ii) We use the angle 6 between the normal to the field 
sheet and the line of sight instead oi (3 — 7r/2 — 6. 

(iii) We make use of GR 8.731 to simplify the expression 
for polarized emissivity. 



The emissivity functions become: 



fi{e) = B^+°'{cose) — Pc^ 



1 -f cos^ ( 
2 cos 6* 



fp{e) = poB'+"(cos( 



\l+aJ - 



1 -f cos^ ( 
2 cos 6* 



(40) 



(41) 



so B^ = Bq (since the field component along the major axis 
is not projected) and rj — cos 9. The emissivity functions are 



For an edge-on field sheet {6 — 7r/2), the argument of the 
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Legendre functions ^ oo and the expressions are replaced ^able 1. Equations for total and polarized emissivity functions 

by: 

for the cases considered in the text. 

fj{n/2) = i3i+«^-V2r(i + |)/r(M^) (42) 
fp = Pofi (43) 

If a = 1, 

// = B^(l + cos^6l)/2 (44) 
fp = poB^sin^6l/2 (45) 



a 




a = l 




r) = 


any rj 


0§ 




0@ 



The total field this time has a two-dimensional Gaussian 



PDF: 

F2uiB)dB = Aexp I I dB 



(46) 



■'0 \ ^^0 , 

We can calculate the emissivities for this Gaussian field dis- 
tribution by integrating B^~^°' over the PDF: 

B^+°'F2D{B)dB = Bl+'^2^r (^^-^^ (47) 

(GR8.310), so: 

ji/fi = jp/fp = 2^r (^) (48) 

If a = 1, this ratio is 2r(2) = 2. The expressions are indeed 
identical to those in equations (bfl) - (p9|). 



3.4 Fields compressed along one axis 

A slightly more corn plex field configuration was studied by 
Hughes et al. (1985) in the context of shocks propagating 
through a jet. We calculate its emissivity in this section in 
order to illustrate the use of the formalism devel oped earlier 
for a three-dimensional case. Hughes et al. (1985) considered 



a cube containing an initially isotropic field which is com- 
pressed to a factor k times its original length along one axis 
(z) and observed in a frame where the plane of compression 
makes an angle e with respect to the line of sight. We take 
the X axis to be in the plane of compression, in which case 
the vector components are: 



ax = 1/fc, bx = 0, cx = 

ay = 0, by = sin e/k, cy ~ cos e 



(49) 



The apparent field position angle xo = 0, so X = X and 
Y = Y. The field components and axial ratio are: 



B% = BoVfc' 

_B|, = _Bo (sin^ e/fc'^ + cos^ e) 

/ • 2 , I 2 2 \l/2 

77 = (sm e + k cos e) ' 



(50) 
(51) 
(52) 



The emissivity functions for a — 1 and the deg ree of 
polarization, p, are given by equations ( p7| ) and ([l8[l, with 
r — 3 (equation E3|): 



/. = |^[2-(l-fc^). 
fp = Po^(l-fc^)co. 



(53) 
(54) 



Po(l 



2-(l-fc2)sin^£ 



(55) 



The last expression is as given by Hughes et al. (1985| ). We 
have also verified that our n umerical results (equ ations - 
|lq) are identical to those of Hughes et al. (1985) for a 7^ 1. 



4 SUMMARY 

We have derived expressions for the synchrotron emissivity 
in Stokes parameters /, Q and U for a three-dimensional, 
disordered, anisotropic field. These are computationally 
straightforward, requiring only two one-dimensional integra- 
tions, which can in turn be used to generate look-up tables. 
In the special case a = 1, the expressions reduce to simple 
analytical formulae. 

The steps in a practical calculation are as follows: 

(i) Start with a cube containing isotropic field, either of 
constant magnitude Bq, or with a Gaussian distribution of 
component rms Bq. 

(ii) Compress and/or shear the cube to get the prescribed 
field configuration (Section B.l; equations p4| - 127 ) . 

(iii) Evaluate the components of the 3 vectors defining 
the edges of the resulting parallelepiped in the observed co- 
ordinate system and normalize them by the volume to give 
the edge vector components ax ,aY ,bx ,by , cx , cy . 

(iv) Work out the apparent field position angle xo from 
equation (|l|). 

(v) Derive the field components along the principal axes, 
and their ratio rj (equations ^ and ^) . 

(vi) Calculate the total and polarized emissivities. The 
relevant pairs of equations for a Gaussian initial field 
strength distribution are summarized in Table 1. 

(vii) For a constant initial field strength, divide by the 
factor r given in equation (^). 

(viii) Rotate by xo to obtain observed Stokes Q and U 
(equations ^ and |2l| ). 

We have shown that the formulae reduce to the 
standard expressions in the literature for one- and two- 
dimensional fields and for an initially isotropic field com- 
pressed along one axis. 
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